Criticality in the 2 + 1-dimensional compact Higgs model and fractionalized insulators 
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We use a novel method of computing the third moment M3 of the action of the 2 + 1-dimensional 
compact Higgs model in the adjoint representation with q = 2 to extract correlation length and 
specific heat exponents v and a without invoking hyperscaling. Finite-size scaling analysis of M3 
yields the ratios (1 + a)/v and 1/u separately. We find that a and u vary along the critical line of 
the theory, which however exhibits a remarkable resilience of Z2 criticality. We propose this novel 
universality class to be that of the quantum phase transition from a Mott-Hubbard insulator to a 
charge-fractionalized insulator in two spatial dimensions. 
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Modelling of strongly correlated systems plays a cen- 
tral role in trying to understand unconventional metallic 
states in cuprate perovskites and other systems, which 
do not conform to the Landau Fermi-liquid paradigm 
JD. One avenue of research attempting to establish a 
theory of non-Fermi liquids in more than one spatial 
dimension, focuses attention on effective gauge theories 
of matter fields representing the charge of doped Mott- 
Hubbard insulators, coupled to compact gauge fields 
emerging from strongconstraints on the dynamics of the 
fermions [| [|, D, g. Compact U(l) gauge fields ex- 
hibit topological defects in the form of monopole con- 
figurations. It has been suggested that the unbind- 
ing of such monopoles may be relevant for spin-charge 
separation in strongly correlated systems ||, M and 
for describing quantum antiferromagnets when fluctua- 
tions around the flux-phase are taken into account 
One often arrives at a description in terms of three- 
dimensional d = 3 compact QED (cQED^). A formu- 
lation of charge-fractionalization in terms of a Z2 lat- 
tice gauge theory coupled to matter fields, has also been 
put forth [ll| . The above provides a link between 

important phenomena in condensed matter physics and 
deep issues in high-energy physics, such as confinement in 
QCD, with which cQED 3 shares two essential features, 
namely confinement and chiral symmetry breaking. 

One lattice model arrived at in this context is the 
compact Higgs model defined by the partition function 
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where N is the number of lattice sites, runs over 
the plaquettes of the lattice, = A^Oj — qAj^, and 
Af+j = SfivxAvAjx- We use the variables (x = + 
1], y = l/[/3+ 1]) in discussing the possible phases of this 



model |12|. In Eq. (Q), r3 is the phase of a scalar matter 
field with unit norm representing holons, A p is a forward 
lattice difference operator in direction /1, while Aj^ is a 
fluctuating gauge field enforcing the onsite constraints 
from strong correlations in the problem. 

When g = 0, the matter field decouples from the gauge 
field. The model has one critical point in the universal- 
ity class of the 3DXY model, y c = 0.688 (y c — 0.75 in 
the Villain- approximation), while the pure gauge theory 
is permanently confined for all values of k [ jl5| . When 
q = 1, Eq. (|l|) is trivial on the line x = 1,0 < y < 1, 
with no phase transition for any value of y. On the line 
< x < l,y = 1 the matter field is absent and the the- 
ory is permanently confined ff^] . The phase-structure for 
q — 2, d — 3 was briefly discussed in Ref. [1^ and subse- 
quently investigated numerically [fl6|| , the phase diagram 
is known, cf. Fig. 5 of Ref. [l6|. When y — > there is an 
Ising transition at x c = 0.5678, when x — > there is a 
3DAY" critical point at y c = 0.688. Moreover, a critical 
line (3 c (k) connects these two critical limits. Above /3 C («) 
the system resides in a decon fined- Higgs phase, while be- 
low Pc(k) it resides in a confined phase. We also note 
that the model in Eq. (|l|) with q = 2 recently was pro- 
posed as an effective theory of a microscopic boson lat- 
tice model exhibiting charge-fractionalized phases [ pd[ . 
The confined phase of Eq. (0) is interpreted as a Mott- 
Hubbard insulating phase, while the deconfined-Higgs 
phase is interpreted as a charge-fractionalized insulating 
phase 

No ordinary second order phase transition takes place 
in the case d = 3,q = 1 (l3[ 0. However, two of us 
have recently shown jOI that when matter is coupled to 
a compact gauge-field in a continuum theory and treat- 
ing the topological defects of the gauge-field in an anal- 
ogous manner to that done in Ref. [l5|, the permanent 
confinement of the pure gauge theory is destroyed. A 
confincment-dcconfincmcnt transition may take place via 
a Kosterlitz-Thouless like unbinding of monopole config- 
urations in three dimensions due to the appearance 
of an anomalous scaling dimension of the gauge-field in- 
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duced by critical matter-field fluctuations The role 
of an anomalous scaling dimension has also been studied 
recently at finite temperature, in pure compact QED in 
d = 3 with no matter fields present | p0[ . In both Refs. 
[E§ £|| the appearance of an anomalous scaling dimen- 
sion is crucial. The authors of Ref. |2(]] recently also 
considered Eq. (Q) with q = 1 numerically [plj, finding 
a recombination of monopoles into dipoles connected by 
matter strings, consistent with Ref. [ ft8| . 

Given the relevance of the case q = 2 to current cen- 
tral issues in condensed matter physics j^, [l0|, O, p"4[ , 



the universality class of the phase transition across the 
critical line for q — 2 warrants attention. We therefore 
compute the critical exponents a and v. Our results i) 
demonstrate that the critical behavior found in the limits 
(3 — > oo (Z2) and the limit k — > 00 (U(l)) are not iso- 
lated points, and ii) on balance suggest that Eq. (|l]) is a 
fixed-line with non-universal a and ^ depending on k) 
rather than exhibiting a Z2- and a XY universality class 
separated by a multicritical point. 

We express Eq. (g) as follows §, fl| 




where S is the Kronecker-delta, D{j — k,m 2 ) = (— A 2 + 
to 2 ) -1 ^jfc, and m 2 = q 2 /3/n. Zo(/3,k) is the partition 
function for massive spin waves and will hereafter be 
omitted. Note the constraint A„ Jj„ = in the func- 
tional integral. Here Qj is the monopole charge on the 
dual lattice site number j, while Jj V are topological cur- 
rents representing segments of either open-ended strings 
terminating on monopoles, or closed loops ]l2] |. 

In the limit (3 — > 00 at fixed k, Eq. (0) takes the form 

{Q}{^} \ j / 

This is the loop-gas representation of the global Z 9 spin 
model in the Villain approximation p2[ . From Eq. ||), 
it is seen that the cases g = 1 and q =/= 1 arc funda- 
mentally different. For q = 1, the summations over 
{Qi} may be performed to produce a unit factor at 
each dual lattice site, eliminating the constraint. Hence, 

/ 2 \ N 

Z= i? 3 (0,e- 27rK ) where i? 3 is an elliptic Jacobi func- 
tion. No phase transition occurs at any value of k for 
q = 1 when [3 — > 00. For q = 2, a phase transition 
survives |r|. In the language of Eq. (||), this crucially 
depends on the presence of the constraint A u Jj V — qQj. 
For g / 1, summing over all values of {Qj} still pro- 
vides a remnant constraint ensuring a theory sustaining 
a phase transition. 

We compute the third moment M3 of the action Eq. 
(|, S = (27rVg 2 )£, Jf v ee (2^/q 2 ) if, with M„ 
given by 

M n = ((H-(H)T). (4) 

Using finite-size scaling (FSS) at the critical point, the 
peaks in M n scale with system size L as L( n ~ 2+a ^ l/ . The 



width between the peaks in M 3 scales as L -1 /", see Fig. 
|l|. Thus both a and v are found from computing one 
quantity without using hyperscaling. Also, FSS of M3 
provides superior quality scaling compared to M2 (spe- 
cific heat), which unfortunately often is marred by sig- 
nificant confluent singularities. We have used multihis- 
togram reweighting p3] | of rawdata and jackknifing for 
error estimation in M3 to perform FSS of the peak height 
and width between peaks. 
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FIG. 1: Generic third moment of action, M3, showing how 
FSS is used to extract a and v. 

Before computing M 3 of Eq. (||), and Eq. (|]), we 
perform benchmark Monte-Carlo simulations (MCS) on 
three well-known models. In Fig. || a) we show FSS 
results for the height of the peaks in M3, defined analo- 
gously to Eq. (|J), for the 3D Ising- and 3Z?XF-models 
for L = 8, 12, 16, 20, 32, 40, 64, with standard Metropolis 
updating. These are limiting cases of Eq. (0) (see be- 
low). Moreover (see below) the 3D Ising spin model is 
dual to the 3D Ising gauge theory (IGT) Q, and we 
have thus also computed M3 for IGT. Any action must 
have a and v identical with those of its dual counterpart, 
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since a and v can be obtained directly from scaling of the 
free energy, and are independent of the degrees of free- 
dom one chooses to describe the system in terms of. Our 
simulations bear this out with precision, cf. Fig. a), 
providing a nontrivial quality check on them. 

The system sizes we have used for the MCS on Eq. 
(|) are L 3 , with L = 8,12,16,24,32,48,64, results are 
shown in Fig. || b) (D-symbols). The allowed MC moves 
using Eq. (|^) are i) insertions of elementary loops made 
of vortex segments Jj V = ±1 and ii) insertions of open- 
ended vortex segments Jj„ = ±q satisfying the constraint 
in Eq. (||). Up to 4- 10 6 sweeps over the lattice have been 
used, with periodic boundary conditions in all directions. 

From Eq. |[ the limit (3 — » oo, n fixed leads to the con- 
straint A^Oj — qAj^ — 2-ivljfj, where Ij^ is integer valued. 
Substituting this into the gauge-field term in Eq. (|l|) , we 
find 
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where = e^xA^ljx 6 Z. For q = 1 the model is 
again seen to be trivial. Since Eqs. (§) and (§) are dual, 
and Eq. (||) is a loop-gas representation of the global 
Z q theory while Eq. (|J) is the Z q lattice gauge theory, 
it follows that the global and local Z q theories are dual 
in d = 3 Hence, the model Eq. (|l]) in the limit 

f3 — ► oo, k fixed, should have a ratio (1 + ot)jv consistent 
with the Z q spin model universality class, if the transition 
is continuous. 

For q = 2, we have performed large-scale MCS and 
FSS analysis of M 3 of the action S in Eq. [l], written as 
S = (3 H$ + k H A , with H x = £[1 - cos(A w -)], cf. Eq. 
(Q). A critical line (3 c (k) separates a confined (/3 < j3 c ) 
and a Higgs-deconfined (/? > /3 C ) state We have 

used L = 8, 12, 16, 20, 24, 32, 40, and up to 9 • 10 6 sweeps 
over the lattice with periodic boundary conditions in all 
directions. The critical line is crossed along the trajec- 
tory (3{n) = f3 c + a (k— k c ), where (/3 C , k c ) is a point on 
the critical line. For the points at which (1 + ot)/v has 
extrema, we use a = (— 1, 1, oo) to check that values for a 
and v are not artifacts of how the critical line is crossed. 

In Fig. b) we show scaling plots of the peaks in M 3 
for Eq. (fy) with q = 2 for various values of k//3 on the 
critical line, Fig. ||c) shows corresponding scaling plots of 
the width between the peaks. From the finite-size scaling 
of the features in M3, Fig. [l], we extract the combination 
(1 + a)/u as well as the exponent 1/u (and hence a) 
along the critical line, the results are shown in Fig. ||. In 
Fig. H c) , we also give values of a obtained directly from 
M3 as well as using (1 + oi)/v together with hyperscaling 
a = 2—di>. We have checked that the extrema in (l+a)/v 
are not changed when the critical line is crossed in three 
very different directions, using a = — 1, a = 1 and a = 00. 

The results seem to rule out that Z2- and XF-critical 
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FIG. 2: Log-log scaling plots of peak-heights in M3. a) 3D 
Ising-, Ising gauge-, and XY-models. b) Eq. (Q) with q = 2, 
= 0.303 (O), «//3 = 1-00 (A), k/0 = 2.88 (a), k/0 = 
3.47 (■), k//3 = 4.05 (x), k/(3 = 5.14 (0). Also shown are 
results for Eq. (0) (□_). c) Scaling plots of width between 
peaks in M3 for Eq. (pi), and Eq. (|l|) for q = 2, legends as in 
b). 



behaviors are isolated points at the extreme ends of the 
critical line. However, from Fig. ||, it is feasible to sug- 
gest two types of universality, Z2 and XY, separated at 
a multicritical point. We believe this to be ruled out 
by the strong deviation in (1 + a)/v from Z2- and XY- 
values at intermediate k//3, which are insensitive to a. 
On balance, we thus conclude that the model Eq. (|l|) de- 
fines a fixed-line theory, rather than exhibiting two scal- 
ing regimes separated by a multicritical point. However, 
the Z2 character of the confinement-deconfincmcnt tran- 
sition persists to surprisingly large values of k/@ on the 
critical line, cf. Fig. 5 of Ref. [l(| Fixed-line theories 
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FIG. 3: a) T = (1 + ck)/k from FSS finite-size of M 3 for 
Eq. (^) for q — 2. b) Same for the exponent u, computed 
directly from M3 (A) and using results for (1 + a)/v with 
hyperscaling (O)- c ) a as computed directly from M3 (A) 
and using hyperscaling (O)- The maximum and minimum in 
a) have been obtained by crossing the critical line along the 
trajectory (3(k) = C + a (re — re c ) with a — 00 (A), a = 1 (■), 
and a = —1 (a) using /3 C = 0.665, k c = 2.125 (max.), and 
/3 C = 0.525, k c = 5.0 (min.). At the maximum, T(a = 00) = 
2.25±0.05, T(o = 1) = 2.23±0.03, T(o = -1) = 2.19±0.06. 
At the minimum, T(o = 00) = 1.30 ± 0.04, T(a = 1) = 
1.31 ±0.05, T(a = -1) = 1.29 ±0.03. Dotted horizontal lines 
indicate ^-values, solid horizontal lines indicate [/(l)-values. 



in 2 + 1 dimensions are known pq |, and non- universal 
exponents imply the existence of marginal operators in 
Eq. ([!]), yet to be identified. 

Recently, Eq. (|l|) with q = 2 was proposed as an ef- 
fective theory for a microscopic model exhibiting a quan- 



tum phase transition from a Mott Hubbard insulator to 
a charge-fractionalized insulator in two spatial dimension 
Q. We thus propose that the zero temperature quan- 
tum phase transition from a Mott-Hubbard insulator to 
a charge-fractionalized insulator |uj jl4| is characterized 
by a fixed-line theory as given in Fig. |p3ut with remark- 
able Z2 resilience. 
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